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We show that, within the seesaw mechanism, an almost decoupled right-handed (RH) neutrino 
species Ndm with mass Mdm ^ 100 GeV can play the role of Dark Matter (DM). The Ndm^s can 
be produced from nonadiabatic conversions of thermalized (source) RH neutrinos with mass Ms 
lower than Mdm- This is possible if a non-renormalizable operator is added to the minimal type I 
seesaw Lagrangian. The observed DM abundance can be reproduced for Mdm 5^^* ~ 10"" Aefi C, 



where Acff is a very high energy new physics scale, 5 
determined by the RH neutrino couplings. 

PACS numbers: 14.60. St, 95.35. +d 



{Md m — Ms)/Md m and ^ < 1 is a parameter 



I. INTRODUCTION 

The results from neutrino oscillation experiments rep- 
resent a success for the seesaw mechanism the sim- 
plest way to understand why neutrinos are massive, yet 
so light compared to all other massive particles in the 
Standard Model (SM). 

Indeed, within the seesaw, the atmospheric and the 
solar neutrino mass scales point to a high energy scale 
~ 10^^ GeV compatible with grand unification and at the 
same time one can understand the observed large mixing 
angles. Moreover, neutrino oscillations support leptogen- 
esis f3], an attractive way to explain the observed baryon 
asymmetry of the Universe and a direct consequence of 
the seesaw mechanism. 

Despite the great progress made in recent years in de- 
riving, especially from leptogenesis interesting con- 
straints on those seesaw parameters that escape the low 
energy experiments investigation, we still lack a way to 
probe the seesaw mechanism. The main obstacle is that, 
for natural choices of the seesaw parameters, the heavy 
right-handed (RH) neutrinos, predicted by the seesaw, are 
not expected to be detected at colliders, because they 
would be either too heavy or too weakly coupled. More- 
over they usually decay very fast disappearing from the 
cosmological lore. If leptogenesis is the right explanation 
of the observed matter-antimatter asymmetry of the Uni- 
verse, produced from the CP violating decays of the RH 
neutrinos, this would be the only reHc trace left over at 
present. 

However, in this paper, we show that a weakly cou- 
pled RH neutrino species can play the role of cold DM. 
The scenario we present differs significantly from the one 
proposed in 0|, where the lightest RH neutrino with a 
O(KeV) mass plays the role of warm DM, and neutrino 
Yukawa couplings are much smaller compared to charged 
leptons and quarks Yukawa couplings. In our model, we 
assume that all RH neutrinos are heavy., with the lightest 
RH neutrino mass not lower than the electroweak scale. 
In this way, the neutrino Yukawa couplings can be of the 
same order as for the other massive fermions. 



II. FAILURE OF THE MINIMAL PICTURE 

The (type I) seesaw mechanism p| is a minimal way 
to explain neutrino masses. The SM Lagrangian is ex- 
tended adding a Yukawa interaction term between three 
RH neutrinos vr and the three left-handed doublets I via 
a Higgs doublet 4> and a Majorana mass term M, 



Y+M 



= lL<t>hvii~]-v'j^MvB. 



(1) 



where h is the matrix of the neutrino Yukawa coupHngs. 

After electroweak symmetry breaking, induced by the 
Higgs VEV V, the Yukawa interaction generates a Dirac 
mass term mo — hv. In the seesaw Hmit, M mo, 
the spectrum of mass eigenstates splits into three light 
neutrinos i^i with masses given by the seesaw formula. 



diag(mi,m2,m3) 



-C/t 



It* 



(2) 



where U is the leptonic mixing matrix, and into three 
heavy neutrinos Ni with masses Mi < M2 < M3. These 
coincide, with very good approximation, with the eigen- 
values of the Majorana mass matrix. 

Neutrino oscillations experiments measure two neu- 
trino mass-squared differences. For normal schemes one 
has to| — 

for inverted schemes one has ; 



^sob 



and ; 



,2 

'atm • 



For mi > matm = yAr 



, -I- Am^ , = 

(0.050 ± 0.001) eV [d| the spectrum is quasidegenerate, 
while for mi < msoi = ^Am^^j = (0.00875±0.00012) eV 
01 it is fully hierarchical (normal or inverted). For defi- 
niteness we will refer to the case of normal schemes but 
all the discussion applies to inverted schemes as well. 

The RH neutrino decays can be conveniently described 
in terms of the decay parameters Ki = TDi/H{T = Mi), 
where Tui are the decay widths. These can be related 
to the neutrino masses introducing the effective neutrino 
masses, defined as fhi = {m}j^mu)u/Mi, such that Ki — 
rhi/rrii,, where m* ~ 1.08 x 10^"^ eV. Assuming Ni to be 
heavier than the Higgs boson. 



2 



from the LEP bound one has Mi > 115 GeV and 
the Ni Hfetimes are then given by 

2 



m^Mf 



5 



TeV 



X 10 



-13 , 



(3) 



Let us now impose that one among the three RH neutri- 
nos species Ni, plays the role of DM particle which we in- 
dicate with Ndm- This implies tom > to — 4 x 10^^ sec, 
where to is the age of the Universe. However, since the 
^DM-decays would produce ordinary neutrinos, a much 
more stringent lower bound comes from neutrino tele- 
scopes Q, 



> a > 1 • 



(4) 



In the range Mdm ~ 10^"^ GeV, the AMANDA Hmits on 
neutrino flux implies a ^ 10^ [9, while in the range 
MjjM ~ 10^~^ GeV, where the atmospheric neutrino flux 
is observed, the lower bound is more relaxed. In any case, 
since strong future improvements are expected from the 
ICE-CUBE experiment, we will leave indicated the de- 
pendence on a in the following discussion |2^. From the 
relation ([3]), this translates into an upper bound on the 
decay parameter Kdm (or equivalently on the effective 
neutrino mass ffiDM) given by 



Kdm (mDM/eY) 



< 



mT 



(5) 



Moreover, imposing that the 7V_DAf-abundance explains 
the measured DM contribution to the energy density 
of the Universe, one flnds a condition on r^M = 
{Nnom /^■y)p'--od, the ratio of the number of Njjm to 
the photon number at the time of the iV£)Af -production, 
occurring at temperatures higher than the electroweak 
phase transition. 



Assuming that the correct value oi tdm is produced by 
some external mechanism, for example from inflaton de- 
cays, a trivial DM model is obtained if the condition 
Eq. (HI is satisfied. Within such a scenario one can in- 
differently identify either A^i or N2 or iVs with Ndm- 
The orthogonal seesaw matrix [l^ . is a useful tool to 
parametrize the Dirac mass matrix niD, such that 



mo = UDll^nDl'^ 



'M 



(7) 



with Dm = diag(mi, m2, TO3) and Dm = 
diag(Mi, M2, M3). The effective neutrino masses 
can then be expressed as linear combinations of the neu- 
trino masses fhi = J2h "^/i \^hi\'^ and one easily obtains 
fhi > mi. Therefore, the upper bound Eq. ^ appHes to 
TOi as well, implying hierarchical light neutrinos. It also 
implies that fl has to be close to the special form 



1 

cos uj sin u! 
— sin u! cos u! 



(8) 



or to those other two obtained by column cyclic permu- 
tation. Therefore, assuming exactly one of these three 
forms for the orthogonal matrix, the condition Eq. ^ 
is fulfilled only assuming some mechanism for the Ndm- 
production based on physics beyond the type I seesaw 
SM extension. Even allowing small deviations from these 
special forms, one undergoes a severe obstacle within the 
type I seesaw. Indeed one can think of different processes 
producing the A^uM-abundance, such as inverse decays 
or scatterings involving the top quark or gauge bosons. 
However, in all cases one has approximately rDM ^ Kdm 
and it would then be impossible to satisfy simultaneously 
the two requirements Eq. ^ and Eq. |[6]). 

Let us consider a particular example that clearly shows 
such a difficulty but that at the same time, as we will see, 
will suggest a solution relying on a simple and reasonable 
extension of the type I seesaw lagrangian. 

We investigate the possibility that the Ndm- 
production is induced by the mixing of Nd m with one of 
the other two RH neutrinos acting as a source, and that 
we indicate with Ns. Notice that Ns has necessarily a 
thermal abundance if the reheat temperature is approxi- 
mately higher than Ms . This is because there cannot be 
more than one RH neutrino species with fhi < m*. 

For definiteness we can assume that Ndm and Ns are 
the two lightest RH neutrinos and hence there are only 
two possibilities: either Mdm = Mi and Ms — M2 or 
vice-versa. In this case A3 does not play any role in the 
A^z?M-production but it is necessary to reproduce cor- 
rectly the neutrino masses. 

This scenario is realized choosing the following form 
for the orthogonal matrix 



— £ 



n' 



(9) 







representing a perturbation, with cos w = 1, of the special 
form in Eq. |[8]). Here the prime index indicates that we 
are reexpressing O into a basis where the RH neutrino 
mass term is still diagonal but in a way that Mdm is 
always the first eigenvalue and Ms the second eigenvalue. 
Notice that we can choose s real and for convenience 
positive. Moreover notice that the choice cosu = 1 is 
not restrictive. Indeed, in any case a value cos uj 7^ 1 
would not be relevant for the DM production but notice 
that it would be important if one simultaneously imposes 
successful leptogenesis from Ns decays, a possibility that 
will be discussed elsewhere 

In order to describe the RH neutrino mixing, it is 
convenient to work in the "Yukawa basis", where the 
Yukawa interaction term is diagonal. This can be diag- 
onalized by mean of a bi-unitary transformation, Dh = 
diag(hA,hB,hc) — VlhUR. The RH neutrino mixing 
matrix Ur can be found considering that it diagonalizes 

h, namely U]^{h'^ h)UR = diag(/i^, /i^). Then, 
from the expression Eq. ([7]), one can see that our choice 
for Vt' simply results into 



3 



cos 6 — sin 9 
sin 9 cos 9 
1 



(10) 



with sin ~ e \/ Ms /Mbm and into 



tlA ^ , riB ^ , ftc — ■ 

V V V 

(11) 

This clearly shows that though does not mix, it 
is necessary to reproduce the atmospheric neutrino mass 
scale. Imposing the condition ([5]), one can see that e has 
to be tiny. Indeed one has 



mi ~ mi + msoi |e|^ 



(12) 



and therefore the upper bound Eq. ^ translates into the 
upper bounds [23| 



mi 
eV 



< 



10 



-33 



TeV 

M 



DM 



< 



10 



-16 



TeV 
M 



DM 



(13) 

This implies a hierarchical light neutrino spectrum and 
a tiny mixing angle between the two lightest RH neu- 
trinos. The description of the production of the N^m- 
abundance proceeds very similarly to the case of light 
active-sterile neutrino oscillations and in particular 
to the case described in [lB|, where transitions occur in 
the non-adiabatic regime as it will prove to be in our case. 
Let us write down the hamiltonian for the two lightest 
RH neutrinos in the Yukawa basis. This will be the sum 
of two terms: a pure kinetic term and a second term ac- 
counting for matter effects described by a potential that 
in the Yukawa basis is diagonal and given by 



Vi 



hj T7(8/fc) 



(I = A,B) 



(14) 



in the approximation of ultrarelativistic neutrinos, imply- 
ing E ^ k and T 3> Ms/3. Notice that in any case for 
T < Ms the A^s-abundance is exponentially suppressed 
and the A''£)M-production would stop anyway. In order to 
further simplify the problem, we also employ a monochro- 
matic approximation where all neutrinos have the same 
mean energy value fc ^ 3T. As usual, we can subtract 
from the hamiltonian a term proportional to the iden- 
tity, irrelevant in neutrino oscillations. Therefore, in the 
Yukawa basis, the relevant hamiltonian can be recast as 



AH 



AM2 
12T 



— cos 29 + {vA 
sin 29 



vb) 



sin 29 
cos 29 — {vA 



VB) 

(15) 

where we defined vi = hj/{4AM^) and AM^ = M|- 
M^ij^j . Approximating cos 26* ~ 1 , one can see that there 
is a resonance at a temperature 



T 




(16) 



only if AM^ < 0, i.e. only if Mi = Ms < Mdm = M2. 
Using the Eq. ifTTj) . Tres can be conveniently recast as 



Tres^ lO^M^Af 




(17) 



If Mdm > 2 Ms one has r^s ^ 10'^ Mdm \/v/Ms- In 
this case, introducing z^cs = Mdm /T^cs — 10~^ a/ Ms/v, 
one can envisage a problem. The Ns's thermalize for 
Zcq — {G/KsY^^ ~ 0.8 [1]. Imposing z^^s > ^cq leads to 
an unacceptably large values of Ms , Md and of the reheat 
temperature. Therefore, unless one assumes an initial 
thermal abundance, one is forced to consider the degen- 
erate limit, for (5 = (Mdm—Ms)/Mdm ^ 1. In this limit 
one now obtains Trcs — lO'^ Mdm S^^^ \/v/Mdm and 
zre. - 10-^ ^Mdm/v. For S < IQ-^^ m^^^TcV, 
this time one can have ^ros ^ ^oq- Therefore, the de- 
generate limit has to be considered as a more attractive 
option. 

Because of the tiny mixing angle the transitions at the 
resonance occur in the nonadiabatic regime. Indeed let 
us calculate the adiabaticity parameter at the resonance. 



1 



29 f 



= sin-^ 29 ■ 



|AM2 



6 T^rcs Htgs 



(18) 



Here H^cs — 1-66 y^Tj^g/Afpi is the value of the ex- 
pansion rate at the resonance. Using the conditions 
Eq. (fT3| and Eq. |[4]), one obtains the upper bound 

7ros < 10-2^ (TeV/McAf)^ • The TV^M-abundance tdm 
can then be calculated as the fraction of iVs's that is 
converted into Ndm ■ This is approximately given by the 
Landau-Zener formula, 



Ndm 
Ns 



(1 







(19) 



Comparing with the condition Eq. ([6]), it is evident that 
neutrino mixing between heavy RH neutrinos cannot pro- 
duce the right Ndm abundance, at least not within a 
minimal type I seesaw extension of the SM. This conclu- 
sion is confirmed by more precise calculations beyond the 
Landau-Zener approximation. 



III. A WAY-OUT FROM 
NONRENORMALIZABLE OPERATORS 

Let us consider the possibility that adding higher di- 
mensional effective operators to the minimal type I see- 
saw Lagrangian Eq. ([l]), while not affecting neutrino 
masses and mixing, enhances the iVDA/-production from 
neutrino mixing. In particular let us consider the follow- 
ing dim-five effective operator 



£effOC^|a>|2iV17VB 



A, 



off 



(20) 
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where $ is the usual Higgs field, A is a dimensionless cou- 
pling matrix and Aos is an unspecified very high energy 
new physics scale that we treat as a free parameter. 

This operator yields a new contribution to 'matter ef- 
fects' into the hamiltonian [l^, that in the Yukawa basis 
can be written as 



rp2 



12 A, 



■Xij. 



(21) 



This result follows from the computation of the temper- 
ature dependent finite real part of the RH neutrino self- 

0: 



energy 



Re[SAr(r)] 



where ni,{q) 



A 



(2^ 



m|)nb(g), (22) 



is the Bose-Einstein distribution 



Sl<!-"l -1 

with u being the four-velocity of the thermal bath. As- 
suming zero Higgs mass one then immediately deduce 
corresponding correction to the Hamiltonian l(2T|) . 

We can reasonably assume that ^ Tics/Aoff. In 
this way in the Yukawa basis the total interaction term 
is approximately still diagonal and with the same eigen- 
values. The relevant hamiltonian describing neutrino os- 
cillations can then be written as 



AH 



off 



AM2 



-vb 



l26I- 



12 T V sin26i + i;^ff 



AB 



VB 



(23) 



where we introduced v'J = A/j/(2 AM^ Aes). No- 
tice that the resonance condition on the temperature, 
Eq. (fT6|l . does not change. However, now the mixing 
angle is different and receives a contribution from the 
off-diagonal terms in H^s, such that sin2^?cff — v^fP . 

Imposing again that mixing is responsible for the DM 
production, since we know that the mixing angle 9 in- 
duced by the Yukawa coupling Ha is by far too small to 
play any role, it can be assumed to be exactly zero. This 
is a good feature since otherwise one could have objected 
that radiative corrections could induce a large value any- 
way, spoiling the stability of Njjm on cosmological scales. 
However, if it is exactly zero, one can invoke some sym- 
metry that protects it from radiative corrections. 

Therefore, the adiabaticity parameter can now be writ- 
ten as 



6 ^rcs -^re 



where we used the Eq. (flG]) for Ties and defined 

1 /4 3 /2 

^ = gj hg /Xab- Using again the Landau-Zener ap- 
proximation for an estimation of the Njjm abundance, 
tndm ~ 7rcs, and imposing again the condition Eq. |(6]), 
we obtain the condition 



(25) 



It is easy to verify that the assumption 
/i^ ^ Tjcs/Acff, translates into a condition 



Ms > 10-2 GeV^y^ (5^/3 /A^j , easily verified ex- 
cept for tiny values of A^^. Notice also that using the 
Eq. in]), one can recast ^ ~ {IO'^/Xab) {Ms/TeYf/*. 
From the condition Eq. I|25p . one then finds in the 
hierarchical case, i.e. Mdm ^ '^Ms, 



Ms 



< 



Aeff y flO^V 

1013 Tevj [XabJ 



TeV, 



(26) 



showing that in order not to satisfy Ms > 100 GeV the 
couplings cannot be too large. On the other hand in the 
more interesting degenerate limit {S <^ 1) one finds 



Mdm > 



A, 



eff 



10^3 TeV 



10- 



Xai 



TeV. 



(27) 



showing, conversely, that in order not to have too large 
values of Mdm the couplings cannot be too small. No- 
tice that too large values log (Mdm /TeV) < 5 8 would 
spoil the cosmologically stability of Ndm , leading to un- 
observed neutrino fiuxes at neutrino telescopes. Indeed 
in this case the nonrenormalizable operator and the mix- 
ing with Ms would induce too fast decays of the Ndm 's 
into Higgs and leptons fl^. For Aes ~ Mqut^ Mpi one 
has then Xab ^ 10-^3^ io_ 'jj^g smallness of Xab can 
be explained in two ways. In the case when Aoff ~ Mqut 
the operator l(20|) can be generated radiatively from the 
coupHng to the GUT scale particles. For example, one 
can assume the Yukawa coupling (with the strength h) 
between RH neutrino, Higgs and heavy (m ^ Mqut) 
fermion. This coupling generates at one loop the op- 
erator l(20|) after heavy fermion is integrated out. The 
values of Xab are, therefore, given by h?'{T^Bs) and, if 
hC^rcs) ^ lO-**"^-^, they come out naturally in the de- 
sired region. Alternatively, if the operator (|20l) is gen- 
erated gravitationally (Aoff ~ Mpi) the smallness of the 
coefficients Xab can be explained in the models where 
the effective value of Mpi in the early universe is dif- 
ferent from its present value (e.g. see Refs.PI). How- 
ever, the consequent decay channels at present should 
be estimated with Xab ^ I. A detailed analysis of the 
constraints from decays will be presented elsewhere 
however, it is remarkable that the mechanism is viable 
for reasonable values of the involved parameters. 



IV. CONCLUSIONS 

We presented a new scenario where the role of DM is 
played by heavy RH neutrinos. The scenario is based on 
a mechanism where the DM RH neutrinos are produced 
through mixing enhanced by the additional presence of 
higher dimensional effective operators into the usual type 
I seesaw Lagrangian. The mechanism rehes crucially on 
the fact that is necessary to convert just a very small 
fraction of the source RH neutrinos into the DM RH neu- 
trinos. In this way the additional operator has the effect 
of enhancing the mixing without spoiling any other suc- 
cessful feature of the type I seesaw mechanism and at 



5 



the same time preserving the DM RH neutrinos stabil- 
ity on cosmological times. A straightforward prediction 
of the mechanism is that the Hghtest neutrino mass has 
to vanish. It also seems quite general that the DM RH 
neutrinos decay and this could lead to signatures in cos- 
mic rays. The recent detected excess of positrons in the 
HEAT and PAMELA experiments have been interpreted 
as due to decaying DM particles with a mass higher than 
300 GeV and a lifetime of approximately tdm ^ 10^^ sec 
[20| . Therefore, our mechanism seems to have the right 
features to explain this excess. These results are quite 
interesting since not only they are fully compatible with 
our model but also because the value for the life time 
corresponds to the saturation of the lower bound Eq. [4] 
from the AMANDA data when Mdm ^ 10^^^ GeV and 



a signal should be expected from the ICE CUBE experi- 
ment. 

It should also be noticed that the special orthogonal 
form Eq. ^ predicted by the mechanism corresponds [2l| 
to a particular sequential dominated model |22||. There- 
fore, the proposed scenario for the solution of the DM 
conundrum restricts remarkably the seesaw parameter 
space, providing a potential smoking gun for the seesaw 
mechanism. 
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